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Abstract 

The build up of polyelectrolyte multilayers (PEMs) was observed by a silicon-on-insulator 
(SOI) based thin film resistor. Differently charged polyelectrolytes adsorbing to the sensor sur- 
face result in defined potential shifts, which decrease with the number of layers deposited. We 
model the response of the device assuming electrostatic screening of polyelectrolyte charges by 
mobile ions within the PEMs. The screening length k -1 inside the PEMs was found to be in- 
creased compared to the value corresponding to the bulk solution. Furthermore the partitioning 
of mobile ions between the bulk phase and the polyelectrolyte film was employed to calculate 
the dielectric constant of the PEMs and the concentration of mobile charges. 
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Introduction. Despite the broad potential applications of polyelectrolyte multilayers (PEMs) a de- 
tailed understanding of the build up process and the resulting basic physical properties is still elusive. 
While the multilayer thickness, the water content, the mechanical properties and the swelling behav- 
ior of different PEMs systems have been extensively studied, their electrostatic properties are still not 
fully determined. PEMs are prepared by the layer-by-layer deposition of polyanions and polycations 
from aqueous solutions. 1 ' 2 During the adsorption process polyanion/polycation complexes are formed 
with the previously adsorbed polyelectrolyte layer 3 leading to a charge reversal. 4 The exchange of 
counterions by the oppositely charged polyelectrolyte could be the reason for the counterion concen- 
tration inside the PEMs to be below the detection limit. 5 Thus, it seems that most of the charges 
within the PEMs are compensated intrinsically by the opposite polymer charges and not by the 
presence of small counterions. Related to the intrinsic charge compensation may be the strong inter- 
digitation between adjacent layers found by neutron reflectometry. 6 ' 7 While the potential of the outer 
PEMs surface is well investigated by electrokinetic measurements, 4 not much is known about the in- 
ternal electrostatic properties like ion distribution and mobility. Using a pH-sensitive fluorescent dye 
the distribution of protons within the PEMs has been determined. 8 Assuming Debye screening and a 
constant mobility for all ions within the PEMs the potential drop within polyelectrolyte films com- 
posed of poly(allylamine hydrochloride) (PAH) and poly(sodium 4-styrenesulfonate) (PSS) has been 
calculated. From these measurements an independent determination of the ionic strength and the 
dielectric constant was not possible. Recent X-ray fluorescence measurements have been promising 
in estimating the ion density profile inside the PEMs giving the total amount of free and condensed 
ions. 9 

Direct measurements of the potential drop inside the PEMs will be best suited for determining elec- 
trostatic properties such as the Debye length or the dielectric constant of the PEMs. The capacitance 
of the PEMs can be measured by electrochemical methods such as AC voltammetry. 10 Another ap- 
proach is the use of field effect devices which allows the determination of the surface potential at the 
sensor /electrolyte interface. Obviously, the surface potential variations measured by such a device 
are strongly dependent on screening effects inside the adjacent phase. The deposition of PAH/PSS 
as well as poly(L-lysine)/DNA multilayers and even DNA hybridization have been detected by such 
devices. 11-13 However, a physical model is needed to relate the quantitative response of the sensors 
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to the dielectric properties and ion mobility inside the PEMs. 

Here we show that a silicon-on-insulator (SOI) based thin film resistor is suited to monitor in real 
time the build up of polyelectrolyte multilayers consisting of the strong polyelectrolyte PSS and the 
weak polyelectrolyte PAH. The sheet resistance of the field effect device is sensitive to variations of 
the potential at the silicon oxide surface. The deposition of the differently charged polyelectrolytes 
results in defined potential shifts, which decrease with the number of layers deposited. Applying a 
capacitor model, the observed decrease can be quantitatively explained by assuming reduced elec- 
trostatic screening by mobile charges inside the PEMs compared to the bulk medium outside. 

Experimental Section. All chemicals including PSS (MW 70,000) and PAH (MW 60,000) were 
purchased from Sigma- Aldrich. Buffers were prepared using ultrapure water (Millipore, France) with 
a resistivity > 18 MQ cm. 5 mg/ml polyelectrolyte solutions were prepared by direct dissolution in 
10 niM Tris buffer at pH 7.5 containing 50 and 500 mM NaCl, respectively. 

The sensor chips were fabricated from commercially available silicon-on-insulator (SOI) wafers (EL- 
TRAN, Canon) using standard lithographic methods and wet chemical etching as described in detail 
elsewhere. 14 The top silicon layer of these wafers was 30 nm thick and slightly doped with boron 
(10 16 cm -3 ). Metal contacts were deposited in an electron beam evaporation chamber (20 nm Ti, 
300 nm Au). After evaporation, the sensor chips were cleaned using acetone and isopropanol. The 
chips were glued into a chip carrier and the contacts were Au-wire bonded to the carrier. Afterwards 
the chips were encapsulated with a silicone rubber to insulate the contacts from the electrolyte solu- 
tion. The sheet resistance of the device is dependent on the potential of the SiO^/PEMs interface 
and was measured as described elsewhere. 15 The potential was then calculated from the sheet 
resistance by a calibration curve of the specific SOI wafer. A flow chamber was mounted on top of 
the sensor and a Ag/AgCl reference electrode was used to control the potential of the electrolyte 
solution and for calibration. The setup and the measurement geometry are shown schematically in 
Figure d First, the sensor was equilibrated in the buffer solution. Next, a calibration measurement 
was performed as shown elsewhere. 15 PAH and PSS solutions were injected twice into the flow cham- 
ber to insure full coverage of the sensor surface, starting with the positively charged PAH. After 
obtaining a stable sensor signal, the chamber was rinsed twice with buffer of the same salt concentra- 
tion as the polyelectrolyte solutions. As soon as a stable signal was obtained, the next polyelectrolyte 
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solution was injected and the procedure was repeated up to 20 times. The sheet resistance of the 
thin film resistor was monitored continuously during the multilayer deposition. In separate experi- 
ments the thickness of the deposited polyelectrolyte films was determined by an ellipsometer (Optrel, 
Multiscope, Berlin, Germany) in electrolyte solution. Care was taken to perform the preparation of 
multilayers as close as possible to the conditions employed for the deposition on the thin film resistors. 
The ellipsometry was carried out in the same buffer solutions used for the sample preparation. 

Results and Discussion. During the layer-by- layer deposition of the polycation PAH and the 
polyanion PSS by alternating buffer exchange, the sheet resistance of the SOI sensor was continuously 
observed. The deposition of the alternately charged polyelectrolytes results in defined responses of 
the sensor (Figure EJ). The adsorption of PAH decreases the sheet resistance corresponding to an 
increased (the potential directly at the SiO^/PEMs interface) as positive charges bind to the 
surface. The subsequent adsorption of the negatively charged PSS increases the resistance and thus 
decreases ^s- As the SOI sensor exhibits also a pH sensitivity, 15 the large potential shift between the 
PAH deposition and the subsequent washing step can be attributed to the pH of the PAH solution 
which is decreased to 6.4 by dissolution of the week polyelectrolyte in the buffer of pH 7.5. The 
deposition of up to 20 monolayers was observable by the field effect device. The potential change 
between adjacent deposition steps A^g was determined from the measured sheet resistance using the 
calibration data and is plotted against the number of adsorbed monolayers as shown in Figure El It 
can clearly be seen that the potential jumps A^s decrease with the number of layers deposited. This 
is in contrast to electrokinetic studies, where the surface potential at the outer PEMs/electrolyte 
interface is measured and the steps remain constant over a large number of deposited layers. 4 The 
observed decrease of A-^g can be explained by adapting a capacitor model. 10,16 The sensor device 
with the adsorbed PEMs is modelled assuming three separate domains: (i) The sensor device, which 
is modelled as a one dimensional silicon/silicon oxide structure characterized by its capacitance per 
area Cg. (ii) The PEMs consisting of N monolayers each of thickness d. (iii) The electrolyte solution 
outside the PEMs, where a diffuse electrical double layer is formed at the PEMs/electrolyte interface 
(Figure 0J). In order to proceed we need to model the charge distribution within the PEMs. In 
general, the layers have a complex charge distribution due to, for example, interdigitation between 
polymers from adjacent layers. This overlapping of layers can even be of the same order as the layer 
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thickness 7 suggesting charge neutralization within the PEMs. 5 Therefore, we consider two limiting 
cases (Figure EJ) : a) separate layers inside the PEMs with uniform volume charge density of p = ±^ 
where a may be regarded as the surface charge density of each layer*, and b) overlapping layers 
such that polyelectrolyte charges are completely neutralized within the PEMs, and thus the only 
uncompensated charges occur at the sensor surface and the PEMs/electrolyte interface with the 
surface charge density of ±|. Electrostatic screening by mobile charges is accounted for using the 
linear Debye-Hiickel (DH) theory both within the PEMs medium (which is characterized by the 
screening length k~ 1 and the dielectric constant e) as well as in the diffuse electrical double layer of 
the electrolyte solution (which is characterized by the bulk screening length Kq 1 and the dielectric 
constant for water e w ) . One can thus identify characteristic Debye capacitances per area of Cp = eeqk, 
and C*d = e w EQKo for the two media, respectively. The two models a and b considered for the limiting 
cases described above yield similar results for the behavior of ^g. On the linear Debye-Hiickel level 
even identical results are obtained, indicating that the precise charge distribution within the PEMs is 
not a critical factor in our model (see Supporting Information). To compare our measurements with 
the model predictions, we calculate A^s — 4>s{N — 1) — ips{N) which can be simplified for nd <C 1 
and even number of layers, N, yielding 

A ^ N>) = (Cs/Cp + C P /C D ) sinh (KNd) + (1 + C S /C D ) cosh ( K Nd) ' ^ 

In Eq. (|TJ) the number of layers, N, appears only in the hyperbolic functions. A value for the screen- 
ing length k _1 inside the PEMs can be obtained from the measured potential change Atps{N) if 
the thickness d of the polyelectrolyte layers at a given ionic strength is known. Therefore the mono- 
layer thickness d was independently measured by ellipsometry. We found an average thickness of 
d = 1.3 ± 0.1 nm for deposition from 50 mM and d = 2.2 ± 0.1 nm for deposition from 500 mM 
bulk electrolyte solution, respectively. Selected ellipsometry data have been crosschecked by X-ray 
synchrotron reflectivity experiments 18 confirming the obtained PEMs thicknesses. As can be seen in 
Figure El the measured Alps can be fitted for nd by Eq. yielding = 6.5 ± 1.0 nm for 500 mM 
bulk solution. A similar fit results in k~ 1 = 6.3 ± 1.0 nm for the build up of the PEMs at 50 mM 



*A mesh size of w 30 nm can be estimated for the adsorbed polyelectrolyte layers. 17 This mesh size is smaller than 
the estimated effective persistence length of PSS 100 nm), suggesting that alternating polyelectrolyte layers form 
with charge densities which are equal in magnitude. 
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bulk solution. Eq. JIJ also provides an estimate for the surface charge density a of the adsorbed 
polyelectrolyte layers, assuming that Cs is much smaller than Cd (1 + Cs/Cd ~ 1). This leads to 
a = 0.020 for the adsorption from 50 mM and a = 0.022 from 500 mM bulk solution. 

The relative dielectric constant e of the PEMs can be calculated from k, given the definition 
k 2 = 2Nxe 2 c/ (esokT) , where N\ is Avogadro's number, e is the elementary charge, c is the ion 
concentration, and k is the Boltzmann constant. For this, one has to determine the concentration c 
of mobile ions inside the PEMs, which is set by the thermodynamic equilibrium between the ions 
in the bulk solution (of concentration Co) and those in the PEMs (see Supporting Information for 
details). It turns out that the dominant factor governing the ion-partitioning in the PEMs/bulk 
electrolyte system is the Born energy change, 
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which arises because of the difference in self-energy of the ions (of radius a) in the PEMs and in 
the bulk leading to the well-known ion-partitioning law c = Coexp(— A/i/fcT). 19,20 Combining the 
preceding relations and the definition of k, the relative dielectric constant e of the PEMs and the 
concentration of mobile ions can be calculated numerically from the values obtained for k. We find 
e = 30 ± 2 and £ — 21 ± 1 for the multilayers adsorbed from 50 mM NaCl and 500 mM NaCl, 
respectively. The corresponding concentration of mobile ions inside the PEMs is estimated to be of 
the order 0.9 ± 0.3 mM and 0.6 ± 0.2 mM, respectively. 

A slightly higher value of £ = 50 ± 10 for the relative dielectric constant of PAH/PSS films has been 
estimated by comparing pyrene fluorescence data of the films with that of various isotropic solvents 
of low molecular weight. 21 Durstock and Rubner 22 found dielectric constants by factor 20 higher for 
PAH/PSS multilayers in water vapor. The large deviation from the values of the present paper is 
not fully understood. A possible explanation could be differences in swelling behavior in water vapor 
and liquid water as shown by neutron reflectometry. 

The different values for £ obtained for PEMs deposited from different salt concentrations can be in- 
terpreted in terms of a different water content of the polyelectrolyte films. A water content of about 
40% was estimated by neutron reflectometry for PAH/PSS films deposited from different salt con- 
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centrations 7 ' 23-25 suggesting that the ionic strength of the solution does not change the water content 
of the PEMs. Assuming an equal water content for both salt concentrations the observed variation 
of the dielectric constant could be ascribed to a different fraction of immobilized to free water within 
the polymer layers as oriented water molecules show a lower dielectric constant. Decreased water 
mobility inside PEMs has already been determined by NMR studies. 26 Comparing measurements at 
different conditions will be necessary to further determine the origin of the dielectric constants of 
PEMs. 

Note that the DH approximation used in the present theoretical model is valid for relatively small 
electrostatic potentials. At room temperature, for symmetrical monovalent electrolytes this yields 
an upper limit of 50 mV, which is typically larger than the potentials measured in our experiments. 
In general a full non-linear Poisson-Boltzmann analysis would be necessary, which, however, is not 
analytically solvable for the present system. The main advantage of the DH approach lies in obtain- 
ing a simple analytical expression for the sensor device functionalized by PEMs enabling a direct 
comparison with the experimental data. 

Conclusion. We were able to show that the recently introduced field effect device based on SOI is 
well suited for the quantitative determination of charge variations at complex interfaces. We apply a 
capacitor model including electrostatic screening by mobile charges within the PEMs to determine 
their dielectric constant as well as the concentration of mobile ions inside the polymer film. The origin 
of the dielectric constants found for PEMs deposited from different salt concentrations will need to 
be addressed further. The presented theoretical description, which is given here for the PEMs, may 
prove useful also for the quantitative analysis of differently functionalized field effect devices. 
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SFB 563 and partially by the French-German Network and by the Fonds der Chemischen Indus- 
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Figure 1: Sketch of the setup and the measurement geometry. Silicon is shown in light grey, silicon 
oxide in dark grey. From top to bottom: native oxide (1-2 nm), conducting top silicon (30 nm), 
buried oxide (200 nm), bulk silicon (675 /zm). A voltage is applied between the source and the drain 
contacts and the resulting current Jsd is measured yielding the sheet resistance of the device. The 
carrier concentration in the top silicon layer is tuned by the backgate voltage U-qg- The potential of 
the electrolyte solution is controlled by a Ag/AgCl reference electrode. A microfluidic device allows 
the rapid exchange of electrolyte solution. 
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Figure 2: The deposition of four PEMs from 50 mM NaCl is shown in detail. The sheet resistance 
was monitored continuously during the build-up process. For each deposition step the polyelectrolyte 
solution was injected twice. When a stable signal was obtained the sensor was rinsed twice with 
buffer. The resistance change AR between adjacent deposition steps of PAH and PSS is indicated 
for two and four adsorbed layers. The potential change between adjacent deposition steps Atp$ wa s 
calculated from AR using the calibration data. The inset displays the subsequent adsorption of 10 
layers. 
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Figure 3: The potential change A^s calculated from the measured change in sheet resistance is 
plotted as a function of the number of adsorbed monolayers N for PEMs deposited from 500 mM 
NaCl. Error bars are determined from the peak to peak noise of the measurement. Applying the 
capacitor model the solid curve was obtained from the fit by Eq. (JTJ). 
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Figure 4: The system is modeled using three separate domains: (i) the Si/Si02 structure, (ii) the 
PEMs, and (iii) the electrolyte solution. The distribution of the immobile charges (the polyelectrolyte 
charges and the silicon oxide surface charge <7i) is shown for two limiting ) separate layers 

inside the PEMs with a volume charge density p and b) overlapping layers inside the PEMs with 
neutralized charges except for the charges at the sensor surface and the PEMs/electrolyte interface 
with the surface charge density of ±|. Volume charges p are displayed in grey whereas surface charges 
a are displayed in black. For case b) the dotted line schematically shows the potential within the 
system. The potential i[)s at the sensor surface is indicated. 
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Theoretical models for the detection of polyelectrolyte multilayers by the 
SOI-based thin film resistor 

Model a: Alternating layers with volume charges 

We model the polyelectrolyte multilayers sensor system as a series of iV alternately positively and 
negatively charged layers with the volume charge distribution p = ±| (Figure 4a, main text). The 
potential in the silicon/silicon oxide structure is assumed to be linear. It is determined by the capac- 
itance per area of the device Cs = §7 which is given by the effective dielectric constant S\ and the 
thickness d\. The charge of the silicon oxide surface is given by U\. In the PEMs, the Debye-Hiickel 
(DH) equation 

d 2 1 

— -7Tll> (x) - K 2 ll) (x) = p (x) 

ax 2 e 

is solved for each single polyelectrolyte layer with p (x) = The screening length inside the PEMs 
is assumed to be k~ 1 and the dielectric constant is e. The diffuse electrical double layer in the 
electrolyte is described by the dielectric constant for water e w and the screening length Kq 1 which is 
equal to the effective double layer thickness. The actual distribution of counterions at the charged 
surface is diffuse and reaches the unperturbed bulk value only at large distances from the surface, 
which is described by an exponential decay of the potential. However, a diffuse double layer behaves 
like a parallel plate capacitor in which the separation between the plates is given by the screening 
length Kq 1 . Thus we can describe the potential within the electrolyte by a plate capacitor with 
the Debye capacitance per area Cd = e w e K . The surface charges a and a L represent the space 
charges within the seminconductor and the electrical double layer of the electrolyte, respectively. The 
potential difference between the bulk semiconductor and the bulk electrolyte is U tot and is set by the 
reference electrode. The potential tp within the N + 2 domains (i), (ii)i, (ii)jv, and (iii) is given 
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by 



(x) = A'x + A, 

V»(ii) n (x) = C n exp (-kx) + D n exp (kx) + ip n (x) , ra = l,...,iV 
^(x) = B'x + B, 

where C n exp( — kx) + D n exp (kx) is the general solution of the homogeneous DH equation and 
ip n (x) is a particular solution of the inhomogeneous DH equation for p (x) = — | (— 1)". We apply 
the following boundary conditions 

V ( i)(0) = 0, V(i)(0) = -^, 

1 

V'(iii) = ^tot, V'(iii) M = — , 



^(n) 1 (xi) - ip(i) (x ± ) = 0, eip' {ii)i (xi) - ex^ji) (xi) = -<n, 
V'(iii) (^jv+i) - ^(ii)^ (»iv+i) = 0, £ w -0' (iii) (xiv+i) - eV(ii )jv (zjv+i) = 0, 

^(ii)„ M - (&n) = 0, £V(ii) n (&n) - ^(ii)^ (&n) =0, u — 2, . . . , N 



with 



X r 



di + (ra-l)d, n = l,...,iV + l, 



xl = d\ + iVrf + k 1 . 



These conditions allow us to calculate the potential — i>(xi) which determines the sheet resistance 
of the device. If we define the Debye capacitance per area Cp = seqk within the polyelectrolyte 
medium we can write as 



4> S (N) 



(Ti + cr eff 



^ sinh (kN<£) + cosh (kN(£) 



U, 



l) N a eS 



(Cs/Cp + C P /C D ) sinh (KNd) + (1 + C S /C D ) cosh (/ciVd) 
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with the effective polyelectrolyte surface charge 

1 — exp (—K>d) 
1 + exp (— nd) 

Model b: Charges overlapping within the PEMs 

The polyelectrolyte multilayers are modeled as overlapping and charges are assumed to be completely 
neutralized within the PEMs. Thus the only uncompensated charges occur at the sensor surface and 
the PEMs/electrolyte interface (Figure 4b, main text). In this case the potential within the three 
domains is given by 

V>(i)(x) = A'x + A, 
V'(ii) ( x ) = C ex P ( — KX ) + D exp (kx) , 
V'(iii) (x) = B'x + B. 

We apply the following boundary conditions 
^i)(0) = 0, V(i)(0) = -^, 

^ 1 

V'(iii) (xl) = Utot, V'(iii) (xl) = —, 

V'(ii) (xi) ~ (xi) = 0, s^Jii) (xi) - £i"0(i) 0&i) = 
V'(iii) (^2) - V'(ii) (^2) = 0, £ w -0(iii) (^2) - e^fii) (^2) 

with 

xi = tZi, 
x 2 = di + iVd, 
xl = di + Nd + K^ 1 . 
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These conditions allow us to calculate the potential ips = ip(%i) an d we can write according to 



^ sinh (uNd) + ^ cosh (nNd) + (u tot - ^(-1)^1 ) 
_l n„ /TCTT^ 7Z KF7\ j_ n _l C S /C D ) cosh (nNd) ' 



t/>s(AH = v — 

} (Cs/Cp + Cp/Cd) sinh (kNJ) + (1 + 



This is the same result which was obtained for model a if a e s in that model is replaced by |. The two 
models represent limiting cases of the charge distribution. Thus in our models the charge distribution 
within the PEMs is not crucial for the values obtained for k^ 1 and e, respectively. For nd <C 1 and 
even numbers of N the potential difference A-^s = ips{N — 1) —ips{N) simplifies to eq 1 (main text). 
The potential ips can now be calculated within the present model as a function of the number N of 
layers deposited as shown in Figure El 
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Figure 5: The potential was calculated within the model for the deposition from a salt concen- 
tration of 500 niM and plotted as a function of the number of adsorbed layers N. The solid lines are 
a guide to the eye. For the calculation the following parameters were assumed: Cs = 2 x 10~ 3 -^j, 
d = 2.2 nm, k^ 1 = 6.5 nm, e = 21, en = —1.25 x 10~ 2 and a = 2.5 x 10~ 2 In the calculation 
the voltage between the bulk semiconductor and the bulk electrolyte was set to 7 mV. 



Ion-partitioning between PEMs and the bulk solution 

One can calculate the concentration of ions inside the PEMs from the thermodynamic equilibrium 
condition with the bulk solution by minimizing the total free energy density (that is free energy per 
unit volume) of the PEMs/bulk electrolyte system as follows. 

The total free energy density consists of the free energy density of the electrical double layer outside 
the PEMs, J-q, and the free energy density of the PEMs medium T?. Within the Debye-Hiickel 

17 



description, the former contribution may be written as Tq = J 7 ^ + jFg orr + (in units of kT) 
comprising the mean-field electrostatic free energy density of the double-layer, J 7 ^, the Debye- 
Hiickel correlation (or excess) free energy density, J 7 ™™, and the self-energy density of the ions, J 7 ^. 
Note that J 7 ™" takes into account the fact that each ion in the electrolyte is surrounded mostly by 
oppositely charged ions, which amounts to the standard correlation free energy expression J 7 ™ 11 = 
— Kq/127t, where k = (8-kI^co) 1 ^ 2 is the bulk inverse screening length with £g = e 2 / (Ane w eokT) being 
the bulk Bjerrum length t. It follows that within the DH approximation, the mean-field free energy 
J 7 ^ is dominated by the entropy of the ions, which is well-approximated by that of an ideal gas of 
particles, i.e. J 7 ^ ~ 2(c lnc — c ) assuming the bulk 1-1 electrolyte concentration of c . Finally the 
self-energy contribution of ions reads ^ elf = 2c £^/2a (in units of kT). 

Similar expressions may be written for the PEMs medium using the inverse screening length k = 
(StiIec) 1 / 2 , the Bjerrum length £b = e 2 / (AneeokT) and the ionic concentration c. In this case, the 
mean-field free energy is again approximated by that of an ideal gas since the net charge within the 
PEMs is assumed to be zero. One thus has T v {c) = 2(clnc - c) - ^/(12tt) + 2c£ B /2a. 

Thermodynamic equilibrium is imposed by minimizing J--p{c) + J^o( c o) with respect to the ion con- 
centration c inside the PEMs, assuming that c + c is constant, which is equivalent to setting equal 
chemical potentials for the two media. One thus finds the ion-partitioning law c = coexp(— A/i/kT), 
where 



The terms proportional to 1/a (with a being the mean ion radius) are due to the self-energy (Born 
energy) difference of the ions in the two media, whereas the terms proportional to the Debye screening 
length come from the DH correlation free energy difference. Since k~ x and Kq 1 are typically much 
larger than a, one can neglect this latter contribution. In fact, an explicit estimate of the dielectric 
constant of the PEMs, e, including the DH correlation free energy shows slightly larger values (up to 
a few percents) as compared with the results reported in the text (obtained based only on the Born 
energy). The difference is however within the experimental error bars. 



tResibois, P. M. V. Electrolyte Theory; Harper & Row: New York, 1968. 
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